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1. Introduction 


The MODIS Airborne Simulator (MAS) is a modified Daedalus Wildfire 
scanning spectrometer which flies on a NASA ER-2 and provides spectral 
information similar to that which will be provided by the Moderate Resolution 
Imaging Spectroradiometer (MODIS), scheduled to be launched on the EOS- AM 
platform in 1998 (King et al. 1992). The principal investigators for the MAS are 
Dr. Michael King (NASA/GSFC, Greenbelt MD), and Dr. Paul Menzel 
(NOAA/NESDIS, Madison WI). 

The Wildfire Spectrometer was delivered to NASA Ames Research Center in 
April 1991. A single visible channel was added and several spectral channels in 
the infrared port were altered to configure the instrument for the FIRE Cirrus-II 
experiment. In January 1992 the modified Wildfire was then further modified to 
become MAS. Beginning in June of 1992, the MAS has been flown in a series of 
experiments that lasted on average 2 to 4 weeks and were held approximately 
every 6 months. 

The MAS spectrometer acquires high spatial resolution imagery in the 
wavelength range 0.55 to 14.3 microns. A total of 50 spectral bands are available 
in this range, and currently the digitizer is configured before each mission to 
record any 12 of these bands during flight. For all pre-1994 MAS missions the 12- 
channel digitizer was configured with four 10-bit channels and seven 8-bit 
channels. The MAS spectrometer is mated to a scanner sub-assembly which 
collects image data with an IFOV of 2.5 mrad, giving a ground resolution of 50 
meters from 20000 meters altitude, and a cross track scan width of 85.92 degrees. 
More details on the MAS sensor and scanner characteristics are shown in Section 
2. A 50-channel digitizer which will record all 50 spectral bands at 12 bit 
resolution is currently under development and is expected to be completed in 
mid 1994. 

In support of the MODIS program, a MAS LEVEL-1B processing system was 
designed and implemented by the MODIS Science Data Support Team (SDST) at 
NASA Goddard Space Flight Center in 1991. The purpose of the processing 
system is to ingest MAS LEVEL-0 aircraft sensor, engineering and navigation 
data, and produce calibrated, geolocated radiances (LEVEL-1B) in a portable 
format (HDF). The system was designed by Liam Gumley, and is currently 
maintained by Paul Hubanks of Research and Data Systems, Corporation. 

The purpose of this document is to describe the characteristics of MAS LEVEL-1B 
data, the calibration and geolocation methods used in processing, the structure 
and format of the LEVEL-1B data files, and methods for accessing the data. 

Questions about MAS LEVEL-1B data should be directed to: 

Paul A. Hubanks, 

MODIS Science Data Support Team, 

Research and Data Systems Corporation, 

(301) 982-3724, phubanks@ltpsun.gsfc.nasa.gov. 
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2. MAS Instrument Summary 


2.1 MAS Specifications 


The MAS spectrometer is mated to a scanning system of the type described by 
Jedlovec et al. (1989). Table 1 summarizes that characteristics of the MAS sensor, 
platform, and scanning system. 

Table 1. MODIS Airborne Simulator Specifications 


Platform 

Altitude 

Ground speed 

Total field of view 

Swath width 

Field of view 

Pixel spatial resolution 

Pixels per scan line 

Scan rate 

Spectral bands 

Spectral range 

Data channels * 

Bits per channel 
Data rate 

Visible calibration 
Infrared calibration 


NASA ER-2 aircraft 
20 kilometers (nominal) 

206 meters per second (nominal) 

85.92 degrees 

37.25 kilometers (at 20 kilometers altitude) 

2.5 milliradians (instantaneous) 

50 meters (at 20 kilometers altitude) 

716 (roll corrected) 

6.25 scan lines per second 
50 

0.55 - 14.3 pm 

12 (selected from 50 spectral bands) 

8 (4 channels @10 bits, 7 channels @ 8 bits) 

246 Megabytes /hour 

Integrating sphere on the ground - none onboard 
Two temperature controlled black bodies onboard 


* Note: MAS data system to be reconfigured to fifty 12-bit channels in mid 1994. 


2.2 MAS Channel Configuration Summary 

In November /December 1991, the modified Wildfire instrument was flown in 
the FIRE Cirrus-II experiment onboard a NASA ER-2 in coordination with other 
aircraft and satellites over the Coffeyville KS field site as well as the TX/LA Gulf 
coast. In January 1992, the modified Wildfire was converted to MAS 
configuration. In June 1992 the MAS was flown over portions of the Atlantic 
Ocean in the region of the Azores during the ASTEX experiment. During early 
1993 the MAS was flown in the southwestern Pacific Ocean during the 
TOGA/COARE and CEPEX experiments. In July 1993 the MAS was flown over 
the northeastern United States during the SCAR-A experiment. During 1994, 
two missions are planned: a Navy ship-tracks experiment (MAST) and a fire 
experiment over California (SCAR-C). 

Although the MAS instrument is a 50 band spectrometer, the data system used 
until mid 1994 has had the capability to record only 12 channels (at 8 bit 
resolution). The tables that follow show the channel configuration of the MAS 
data system for missions flown through July 1993. The recording of particular 
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bands in specified channels is arbitrary and is determined before each experiment 
by the principal investigators (scientists) for each mission. The X column shows 
the wavelength at 100% response, and the AX column shows the bandwidth at 
50% response of the MAS band recorded in each data channel. Note that channel 
1 is currently used as a "bit bucket" to store the two 'extra' bits for channels 9-12. 
In cases where laboratory measurements of the spectral response for a band were 
not available, Gaussian curves were fit to the predicted central and 50% response 
points. Plots of MAS spectral responses for each experiment are in Appendix A. 

Table 2(a,b). MODIS Airborne Simulator Channel Configuration 


(a) 


Channel 

FIRE (10/91) 
X (n.m) AX 

ASTEX (6/92) 
X (|im) AX 

TOGA/COARE 

(1/93) 

X (|i.m) AX 

Number 
of bits 

1 

2 

0.681 

0.009 

0.664 

0.053 

0.664 

0.053 

8 

3 

1.617 

0.056 

0.875 

0.041 

0.875 

0.041 

8 

4 

1.933 

0.048 

0.945 

0.043 

1.621 

0.057 

8 

5 

2.088 

0.047 

1.621 

0.057 

1.830 

0.050 

8 

6 

2.139 

0.047 

2.142 

0.047 

2.142 

0.047 

8 

7 

3.748 

0.144 

3.725 

0.151 

3.725 

0.151 

8 

8 

4.695 

0.144 

13.952 

0.517 

13.952 

0.517 

8 

9 

4.539 

0.144 

8.563 

0.395 

8.563 

0.395 

10 

10 

8.800 

0.339 

11.002 

0.448 

11.002 

0.448 

10 

11 

10.950 

0.500 

13.186 

0.351 

13.186 

0.351 

10 

12 

11.950 

0.500 

12.032 

0.446 

12.032 

0.446 

10 




(b) 





Channel 

CEPEX (3/93) 

SCAR-A (7/93) 

Number 








of bits 




X (tarn) 

AX 

X Oim) 

AX 




1 

2 

0.664 

0.055 

0.547 

0.043 

8 



3 

0.875 

0.041 

0.664 

0.055 

8 



4 

1.623 

0.057 

0.875 

0.041 

8 



5 

1.880 

0.050 

0.945 

0.043 

8 



6 

2.142 

0.047 

1.880 

0.050 

8 



7 

3.725 

0.151 

2.142 

0.047 

8 



8 

13.952 

0.517 

3.725 

0.151 

8 



9 

8.563 

0.396 

8.563 

0.395 

10 



10 

11.002 

0.448 

11.002 

0.448 

10 



11 

13.186 

0.352 

13.186 

0.352 

10 



12 

12.032 

0.447 

12.032 

0.447 

10 



2.3 MAS Experiment Acronyms 

FIRE First ISCCP Radiation Experiment (10/91 - 12/91) 
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ASTEX 

TOGA/COARE 

CEPEX 

SCAR-A 

2.4 Glossary 

Ames 

Band 

Blackbody 

Browse 

Calibration 

Channel 

Data system 

ER-2 

Exabyte 

Exabyte format 
Flight track (line) 
FTP 

Geolocation 

GIF 

HDF format 

IFOV 

INS 

Intermediate 

IR 

ISCCP 

LEVEL-0 

LEVEL-1A 

LEVEL-1 B 

MAS 

Metadata 

MODIS 

Navigation data 

NIR 

PCM 

Scanner 

SGI 

Spectrometer 

VIS 


Atlantic Stratocumulus Transition Experiment (5/92 - 6/92) 
Tropical Ocean Global Atmosphere /Coupled 
Ocean- Atmosphere Response Experiment (1/93 - 3/93) 
Central Equatorial Pacific Experiment (3/93 - 4/93) 

Sulfates, Clouds and Radiation - America (7/93) 


NASA Ames Research Center 
One detector /filter assembly 

Radiation source, emittance=absorptance=l, reflectance=0 
Uncalibrated reduced resolution image product 
Conversion from sensor digital counts to physical units 
Area in LEVEL-0 format reserved for data from 1 band 
MAS digital flight data digitizer /recorder onboard ER-2 
NASA high altitude research version of U-2 spy plane 
8 mm tape drive for digital data storage 
Storage format for ER-2 Exabyte data system 
Portion of aircraft track where heading is constant 
File Transfer Protocol (Internet) 

Computation of earth coordinates of a sensor measurement 

Graphics Interchange Format 

Hierarchical Data Format 

Instantaneous Field Of View 

Inertial Navigation System onboard ER-2 

(format): Data format created by unpacking LEVEL-0 data 

Infrared: Wavelengths of 3.7 microns and longer 

International Satellite Cloud Climatology Project 

Unformatted, raw, packed, unprocessed sensor data 

Unpacked and reformatted sensor data 

Calibrated and geolocated sensor data 

MODIS Airborne Simulator 

Descriptive information pertaining to a dataset 

Moderate Resolution Imaging Spectroradiometer 

Record of the position and orientation of a sensor platform 

Near-Infrared: Wavelengths from 0.7 to 3.7 microns 

Tape drive used for old ER-2 data system 

MAS scan mirror and blackbody target assembly 

Silicon Graphics Workstation 

MAS optics and detector assembly 

Visible: Wavelengths from 0.4 to 0.7 microns 
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3. MAS Data Calibration 

3.1 Visible/Near-Infrared Band Calibration 


Calibration coefficients for the MAS visible (X < 0.7 |im) and near-infrared (0.7 
(im < X < 3.7 pm) bands are derived from integrating spheres used pre and post 
launch. The transformation from digital count to radiance is defined by 


Radiance 


where 


( Count - ( Gain • CBB aV g ) ) • Slope 


( 1 ) 


Radiance 

Count 

Gain 

CBB a vg 

Slope 


= computed radiance, 

= recorded count, 

= channel gain, 

= running average of 30 previous cold black body counts, 
= slope coefficient, determined for each band. 


The radiance units are Watts per square meter per steradian per micron (W m -2 
sr _1 pm -1 ). Integrating sphere calibrations may be carried out several times 
during a MAS flight experiment. There may be several sets of calibration slopes, 
each applied to data from a discrete time period, during a single experiment. For 
more details, see Arnold et al. (1994a, 1994b). 


3.2 Temperature Correction of Visible/Near-Infrared Bands 


MAS testing in a cold chamber at NASA Ames Research Center has shown that 
some MAS VIS/NIR bands are subject to sensitivity changes as a function of 
temperature. This is important, since at cruise altitude (20000 meters) the MAS 
instrument environmental temperature is typically around -35°C. The cold 
chamber data is used to develop equations which predict the sensitivity of the 
VIS/NIR channels at given temperatures. Temperature sensors mounted on the 
MAS in-flight are then used as a reference for computing the sensitivity 
correction. For example in ASTEX data it was found that only channels 5 and 6 
were sensitive to this phenomenon. For these two channels Count=C e ff, for all 
other VIS/NIR channels Count=C rec . 

The form of the equation for C e ff during ASTEX was: 

Count = C e ff = Gain • C c bb + ( C rec - ( Gain • C c bb ) ) / ( a] T mas + b\ ) (2) 

where 

C e ff = temperature corrected effective digital count. 

Gain = channel gain, 

C c bb = recorded in-flight digital count for cold black body. 
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C r ec = recorded original in-flight digital count, 

a l/ b l = coefficients determined for each band, 

T ma s = temperature of cold black body (degrees Celsius). 

The coefficients ai and bi are determined for each temperature sensitive channel 
in each MAS experiment. This correction is applied before the linear 
transformation to radiance units. For more details on the corrections used in 
various field experiments, see Arnold et al. (1994a, 1994b). 

3.3 Infrared Band Calibration 

Calibration data for the MAS infrared (IR) bands (X > 3.7 pm) are obtained during 
flight from two blackbody sources which are viewed during every mirror scan. 
The first (cool) blackbody is usually left to float at ambient temperature, 
although it can be temperature controlled. The second (warm) blackbody is 
maintained at a higher temperature than the cool blackbody, and the 
temperature level is set depending on the desired scene temperature range. 
Thermistors measure the temperature of the blackbodies during flight. The 
equivalent Planck radiances for the blackbodies (assuming unit emissivity) are 
computed using a sensor weighted integral of the Planck function of the form. 

R(\,T) = 1 B(X,T) S(X) dX (3) 

J S(X) dX 

over the wavelength range of the spectral band where 

R(A/T) = sensor weighted equivalent Planck radiance, 

B(X,T) = Planck radiance at wavelength and blackbody temperature T, 

S(X) = sensor spectral response. 

The Planck radiance in Watts per square meter per steradian per micron (W nr 2 
sr _1 pm -1 ) is computed by the equation 


B(X,T) = 

io- 6 . Cj 


X 5 . ( exp (C 2 / ( X . T ) ) - 1) 

where 

Cl = 

2 . h . c 2 = 1.1910439 . 10‘ 16 W nr 2 

c 2 = 

( h . c ) / k = 1.4387686 . 10‘ 2 m K, 

X 

wavelength in meters. 

T 

temperature in degrees Kelvin. 


Once the equivalent Planck radiances for the blackbodies are known, the 
calibration slope and intercept are computed by the relationships 
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Slope 


r (^/T)bb2 " R(X,T) B bi 

Count BB2 " Count BB1 


(5) 


Intercept = ( R(X,T) BB1 . Count BB2 ) - ( R(X,T) BB2 . Count BB1 ) (6) 

Count BB2 - Count BB1 

where 

R(X,T) BB i and R(X,T) BB2 are the equivalent Planck radiances for the cool 
and warm blackbodies respectively, 

Count BB i and Count BB2 are the digital radiance counts for the cool and 
warm blackbodies respectively. 

The calibration from digital counts to radiance is then a linear transformation of 
the form 

Radiance = Slope • Count + Intercept (7) 

The radiance units are Watts per square meter per steradian per micron (W m -2 
sr -1 |j.m _1 ). This procedure is performed for every IR channel on every scanline. 
No averaging of MAS blackbody data is done. 

Conversion from IR radiance to Planck equivalent temperature or 'brightness 
temperature' may be done by inverting the Planck equation. The inverse 
equation is of the form 

T(X,B) = C 2 (8) 

X.log e (C 1 /(x5.B(X / T).106) + 1) 

where 

T(X,B) = brightness temperature in degrees Kelvin, 

C 2 = ( h . c ) / k = 1.4387686 . 10" 2 m K, 

Q = 2 . h . c 2 = 1.1910439 . 10' 16 W nr 2 , 

X = wavelength in meters, 

B(X,T) = Planck radiance in W nr 2 sr -1 pm -1 . 
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4. MAS Data Geolocation 

Navigation data onboard the ER-2 is recorded by an Inertial Navigation System 
(INS). This system has a dedicated tape recorder, which is separate from the 
MAS tape recorder. Before February 1992, INS data was obtained from the INS 
recorder itself. This data was unpacked and formatted by staff at Ames Research 
Center, and distributed as a separate file along with the MAS data. Since 
February 1992, the new MAS Exabyte data recorder has had the capability to 
record data from the INS system at the same time as it records MAS data, thus 
giving one integrated dataset which contains both MAS and navigation 
information. In the old INS recording system, aircraft position and attitude were 
updated every 5 seconds. In the new recording system, the update interval is 1 
second. 

The geolocation algorithm is only applied to those portions of a flight where the 
aircraft flew straight, level, and on a constant heading. Data recorded during 
ascent, turns, and descent is not processed to LEVEL-1B. Ascent data is not 
processed since the scanner is not looking at nadir, and thus the imagery is 
distorted. During turns, the scanner roll correction system cannot compensate 
for high roll angles, and the imagery is distorted so that black bands appear on 
either side of the flight track. During descent, the scan mirror optics become 
covered with condensation, and the image information is lost. 

The straight line flight track start and end times are identified by an automated 
program which analyzes altitude, heading, pitch, and roll. Once the start and 
end times for all flight tracks are determined, linear regressions for aircraft 
latitude, longitude, heading and altitude versus time are computed. This is the 
simplest way to represent the INS data set when it is referenced to the MAS data 
set, which is sampled at a higher rate (6.25 scanlines /second). To geolocate a 
given MAS flight track, the MAS start time and scanline number at the 
beginning of the flight track are determined. These are used as a reference for the 
rest of the flight track, since the time code recorded by the MAS is stored as 
integer seconds. The scanline number and scan rate are used to determine the 
time at subsequent scanlines in the flight track by 

t = t 0 + ( S - S 0 ) / r (9) 

where 

t 0 = time of the first scan line in a flight track, 

S = scan line number of a data line in the flight track, 

S 0 = scan line number of the first data line in the flight track, 

r = scan rate (6.25 scans/ second). 

Once time for a given scanline is computed, the linear regression relationships 
are used to compute aircraft latitude, longitude, heading and altitude at that 
time. The scan angle is determined from the pixel number by the relationship 
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DD = 

K . (E-E n ) 

(10) 

where 



II II II n 

O - 

O M w m 

scan angle (displacement from aircraft nadir point, degrees), 
angular distance between IFOVs (85.92 degrees / 715 steps ), 
pixel number, 

pixel number at nadir (358.5). 

The distance from the aircraft nadir point to the pixel is then given by 


D 

Z . tan ( DD ) 

(11) 

where 



Z 

aircraft altitude. 



Once the distance and azimuth from nadir to the pixel are known (azimuth is 
computed from the heading), trigonometric transformations are used to 
compute the latitude and longitude for every 10th pixel on a scanline (pixels 1, 
10, 20, 30, 690, 700, 710, 716). The reason for geolocating every 10th pixel on a 

scan line is to save space in the output data set. Solar zenith and azimuth angles 
and sensor zenith and azimuth angles are also computed for every 10th pixel. 
More details on the MAS geolocation algorithm are given in Jedlovec et al. 
(1989). 
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5. Obtaining MAS LEVEL-1B Data 
5.1 Distribution Contact Points 

MAS LEVEL-1B data are currently distributed by the MODIS Science Data 
Support Team (SDST). However it is expected that by the end of 1994 the NASA 
Langley and Goddard DA AC’s (Distributed Active Archive Center) will be the 
distribution points for MAS data. At that time users will need to contact the 
appropriate DAAC for specific media and format information. MAS LEVEL-1B 
data from the FIRE II Cirrus and ASTEX missions will be distributed by the 
Langley DAAC. MAS LEVEL- IB TOGA/COARE data will be distributed by the 
Goddard DAAC. Distribution points for other experiments have yet to be 
determined. 

MODIS Science Data Support Team 

Paul A. Hubanks 

Research and Data Systems Corporation 
7855 Walker Drive, Suite 460 
Greenbelt MD 20770 
Phone: (301) 982-3724 

Internet: hubanks@ltpsun.gsfc.nasa.gov. 

Langley DAAC 

User and Data Services 
Mail Stop 157-B 

NASA Langley Research Center 
Hampton VA 23681 
Phone: (804) 864-8656 

Internet: userserv@eosdis.larc.nasa.gov 

Goddard DAAC 

DAAC User Support Office 
Global Change Data Center 
Code 902.2 

NASA Goddard Space Flight Center 
Greenbelt MD 20771 
Phone: (301) 286-3209 
Internet: daacuso@daac.gsfc.nasa.gov 

5.2 SDST Distribution Tape Format 

MAS LEVEL-1B data distributed by the MODIS SDST are stored on Exabyte 8500 
8mm tape. The tapes are formatted as simply as possible, so they may be read on 
any other system without the problems imposed by system dependent header 
files, block lengths, and formatting. Each tape contains all the flight lines for one 
MAS flight (one day). The number of flight lines per flight varies, but is 
generally between 10 and 20. The volume of data varies, but is generally 1 to 3 
gigabytes per flight. 
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The archive tapes are created by writing each output data file (1 straight-line 
flight track) to tape in fixed-length blocks of 16384 bytes, in time ascending order. 
One end-of-file (EOF) mark is written at the end of the data blocks for each file, 
and an extra EOF is written at the end of the data on the tape. The last block of 
each file has good data at the start of the block and unused bytes (filled with null 
characters) at the end. Information on the length of the file is encoded in the 
header when the file is created. No file name, protection, or ownership 
information is written onto the archive tape. All information necessary to 
identify the file is stored in the file itself. 

To read the data files from tape, system utilities should be used. For example on 
a UNIX system, the following command could be used to read the first file on the 
tape: 

% dd if =/dev/mt/tps0d3nrnsv. 8500 of =outf ileOl . hdf ibs=16384 

Note that "/dev/mt/tps0d3nrnsv. 8500" is the tape device driver. This is system 
dependent. However, the user should ensure that the device driver is for an 
Exabyte 8500 (not 8200) drive, and handles fixed length blocks which are 16384 
bytes long. The output file name is a sample only: the user could specify any 
output file name. 

To read the first straight-line flight track file on a DEC VAX/VMS system, the 
following commands could be used (note that any VMS utility which reads tapes 
with fixed-length blocks is suitable): 

$ allocate mual : tape 

$ mount/foreign/over=own/blocksize=16384/recordsize=16384 tape 
$ set magtape /rewind tape 
$ copy tape: outf ileOl . hdf 
$ dismount tape 
$ dealloc tape 


Tape Format Summary : Fixed length blocks , 16384 bytes per block , one EOF after each 
file , two EOFS at end of tape. Each file represents one flight line. 
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6. How to decode MAS Level-IB Data 
6.1 Introduction to HDF 

MAS LEVEL- IB data are stored in Hierarchical Data Format (HDF). The data are 
organized into individual straight-line flight tracks. Each MAS LEVEL-IB HDF 
file contains one straight-line flight track, and defines a granule of data. In 
general there are 10 to 20 flights (or flight days) in a given experiment and 10 to 
20 flight tracks in a given flight. 

HDF is developed and maintained by the National Center for Supercomputing 
Applications (NCSA), and has been selected as the standard data format for all 
Earth Observing System (EOS) data products. HDF is a multi-object file format 
for the transfer of graphical and numerical data between machines, and allows 
the user to create, access, and share scientific data in a form that is self-describing 
and network-transparent. "Self-describing" means that a file includes 
information defining the data it contains. "Network-transparent" means that a 
file is represented in a form that can be accessed by computers with different ways 
of storing integers, characters, and floating-point numbers. HDF files also allow 
direct access, so that a small subset of a large dataset may be accessed efficiently, 
without first reading through all the preceding data. It should be noted that HDF 
files should only be accessed through the HDF library of subroutine and function calls 
(from FORTRAN or C). HDF is not a data format which can be unpacked easily 
by knowing word locations, byte ordering and so on. With HDF, the user is 
insulated from these details, so that differences in system specific storage details 
are transparent. 

MAS LEVEL-IB data are stored in a subset of HDF known as HDF/netCDF. 
Initial versions of the MAS LEVEL-IB data were written in netCDF (Network 
Common Data Form) format, developed by the University Corporation for 
Atmospheric Research (UCAR). In HDF Version 3.3, the netCDF interface was 
integrated with the HDF software library. All MAS LEVEL-IB data is now 
written in HDF/netCDF to comply with the EOS data format standard. 

From the user's perspective, this means that while the data is stored in 
HDF/netCDF, the data may be extracted using netCDF-like software calls, from 
either FORTRAN or C. This is the preferred means of extracting data from MAS 
LEVEL-IB data files. While it is possible to use HDF software calls alone, this 
method is not described in this document. 


Data Format Summary: MAS data is stored in HDF/netCDF , and may be extracted 
using netCDF-like ca lls to the HDF library. 
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6.2 Obtaining HDF libraries 


In order to extract data from MAS LEVEL-1B data files, the HDF software library 
must be compiled and installed. The source code is available from NCSA by 
anonymous FTP on Internet. The following commands may be used on a Unix 
system to retrieve the source code: 

% ftp ftp.ncsa.uiuc.edu 
Name : anonymous 
Password: guest 
ftp> cd HDF/HDF3 . 3r3 / tar 
ftp> binary 

ftp> get HDF3 . 3r3 . tar . Z 
ftp> quit 

% uncompress HDF3 . 3r3 . tar . Z 
% tar xvf HDF3 . 3r3 . tar . Z 

The instructions provided in the README files should then be followed to 
compile the source code and construct the library. If problems are encountered, 
please contact HDF user support by electronic mail at hdfhelp@ncsa.uiuc.edu. 

Table 3. Platforms on which HDF 3.3 Release 3 has been tested by NCSA 


Platform 

HDF base library 

HDF/netCDF 

Convex/ConvexOS 

X 

X 

Cray Y-MP/UNICOS 

X 

X 

DEC Alpha /OSF 

X 

X 

DecStation / MIPSEL 

X 

X 

Fujitsu VP/UXPM 

X 

- 

HP/UX 9.01 

X 

X 

IBM PC - MSDOS 

* * * 

* * * * 

IBM PC - Windows 3.1 

* * * 

* * * * 

IBM PC - Windows NT 

* * * 

* * * * 

IBM/RS6000 

X 

X 

Intel i860 

X 

- 

Mac/MacOS 

X 

* * 

NeXT/NeXTSTEP 

X 

X 

SGI/IRIX4 

X 

X 

Sun4/Solaris 

X 

X 

Sun4/ SunOs 

X 

X 

VAX/VMS 

X 

X 


X Successfully tested. 

* * The netCDF interface works but old netCDF XDR files cannot be read. In 

addition, the FORTRAN interface and ncdump utility have yet to be 
finished. 

*** There is no FORTRAN support for either PC version of HDF 3.3r3. 

* * * * The netCDF part of the HDF / netCDF merger is not working correctly. 
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6.3 Decoding and Extracting MAS Level-IB Data 

To make it easier to get started using the HDF library, a compiled version for a 
Silicon Graphics (SGI) platform is available via anonymous FTP at 
ltpiris2.gsfc.nasa.gov. In the future, compiled libraries for several different 
platforms will be available. The user should log in with username anonymous 
and password guest. The files are stored in the directory 

pub/MAS /hdf _ne tcdf / SGI__IRIX4 . 05 . 

The following files are provided (A = ASCII mode, B = binary mode in FTP): 


libdf.a 
libnetcdf.a 
ncdump 
netcdf.h 
netcdf.inc 
simple 
simple. f 
simple. sh 


compiled HDF software object library (B) 
compiled HDF/netCDF software object library (B) 
executable utility program to dump HDF files as formatted text (B) 
C header file needed to call the HDF/netCDF library (A) 

FORTRAN include file needed to call the HDF/netCDF library (A) 
executable utility program to read MAS LEVEL-IB HDF files (B) 
FORTRAN program to read MAS LEVEL-IB HDF data (A) 
script to compile the FORTRAN program simple.f (A) 


These files provide the tools to begin extracting data from the MAS LEVEL-IB 
HDF data files using the FORTRAN or C netCDF interface. For further 
information on programming using the netCDF FORTRAN or C interface 
consult the netCDF User's Guide, available in PostScript form via anonymous FTP 
from UCAR as shown: 


% ftp unidata.ucar.edu 
Name : anonymous 
Password: guest 
ftp> cd pub/netcdf 
ftp> binary 
ftp> get guide. ps . Z 
ftp> quit 

% uncompress guide. ps . Z 

Note that if you downloaded the HDF3.3r3 source code as previously described, a 
copy of the netCDF User’s Guide is included in the distribution. 

The following FORTRAN program (simple.f described above) gives a simple 
demonstration of how to extract MAS radiance data from a LEVEL-IB HDF file 
using the netCDF interface. While it does not exercise many of the features of 
the netCDF interface, it does demonstrate the basic functions of opening the file, 
defining the variables to be read, and extracting a portion of the selected data. 
Please note that the netCDF User's Guide is indispensable when developing code to 
read and/or write HDF/netCDF files. 
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c a simple program to demonstrate how to read a MAS HDF file 

c Liam Gumley, RDC, 03-SEP-1992 

c include the netcdf definitions 

include ' netcdf . inc ' 

c set up necessary data types for netCDF (note that you should 

c leave these as the default type for your compiler) 

integer cdfid, rcode, dataid, start ( 3 ), count ( 3 ) 

c set up types for variables and attributes 

integer*2 value 
real*4 scale ( 12 ) 
character string*72 
character name* 7 2 

c set netCDF error options 

call ncpopt ( NCVERBOS + NC FATAL ) 

c open the netcdf file 

write (*,'(' ' Enter HDF file name : ''$)') 

read ( * , ' ( a7 2 ) ' ) name 

cdfid = ncopn ( name , NCNOWRIT, rcode ) 

c get the variable id for the desired variable 

dataid = ncvid( cdfid, ’ CalibratedData ' , rcode ) 

c get the scale factor values (from attribute scale_factor ) 

call ncagt ( cdfid, dataid, ' scale_factor ' , scale, rcode ) 

c get the units text description (from attribute units) 

call ncagtc ( cdfid, dataid, 'units', string, 72, rcode ) 

c set the pixel, channel and record counters 

count ( 1 ) = 1 

count ( 2 ) = 1 

count ( 3 ) = 1 

c get the start record, channel and pixel coordinates 

20 write(*,*) 'enter record, channel, pixel' 

read ( * , * ) start ( 3 ), start ( 2 ), start ( 1 ) 

c get the hyperslab of data (one value in this case) 

call ncvgt ( cdfid, dataid, start, count, value, rcode ) 

c write the resulting radiance, after rescaling to a real number 

write (*,40) real ( value ) * scale ( start ( 2 ) ), string ( 1 :29 ) 

40 forma t ( ' Radiance = ’, fl2.3, lx, a29 ) 

go to 20 
end 
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7. Structure and Contents of MAS LEVEL-1B Data 
7.1 Data Structure Summary 


Each individual MAS LEVEL-1 B file contains calibrated, geolocated radiances for 
all MAS channels for one straight line flight track. Each file also contains 
various ancillary data such as instrument configuration, and time and date 
codes. When the MAS LEVEL-1 B data structure was designed, a major concern 
was that all LEVEL-0 engineering and ancillary input data be maintained in the 
output data set. For this purpose, all LEVEL-0 parameters (except the raw digital 
image counts) are retained in the LEVEL-1B structure. In this way, the health 
and performance of the sensor may be monitored, and the data may be reversed 
to LEVEL-0 if necessary. 

The MAS LEVEL- IB data structure was designed to be easy to understand and to 
reflect the physical characteristics of the data. All variables, dimensions, 
attributes, and data are named descriptively to aid interpretation. The following 
listing was created from an actual MAS LEVEL-1B HDF file using NCDUMP (see 
Section 7.2). Due to size constraints only the dimensions, variables, and global 
attributes are listed here. The data portion is not shown. 

netcdf 92169-08 { 
dimensions : 

Time = UNLIMITED ; // (5367 currently) 

NumberOf Channels =12 ; 

NumberOf Pixels = 716 ; 

HeaderLength = 6160 ; 

AnchorlndexSize = 73 ; 

variables : 

float Lef t50%ResponseWavelength ( Number Of Channels ) ; 

Left50%ResponseWavelength: units = "microns" ; 
float CentrallOO%ResponseWavelength (NumberOfChannels) ; 

CentrallOO%ResponseWavelength: units = "microns" ; 
float Right 50 %ResponseWavelength (NumberOfChannels) ; 

Right50%ResponseWavelength: units = "microns" ; 
float SolarSpectrallrradiance (NumberOfChannels) ; 

SolarSpectrallrradiance : units = "watts/meter2 /micron" ; 
SolarSpectrallrradiance : title = "Sensor Weighted and Orbit 
Corrected Solar Spectral Irradiance" ; 

short AnchorPointlndex (AnchorlndexSize) ; 

char DataSetHeader (HeaderLength) ; 

short DataFrameStatus (Time) ; 

short RunNumber (Time) ; 

long ScanLineCounter (Time) ; 

long Thumbwheel Switches (Time) ; 

short ScanRate (Time) ; 

ScanRate : units = "hertz" ; 

ScanRate : scale_f actor = O.lf ; 
long GreenwichMeanTime (Time ) ; 

long MasTime (Time) ; 
long InsTime (Time) ; 
short S-Bendlndicator (Time) ; 
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short Aircraf tRollCount (Time) ; 
long YearMonthDay (Time) ; 
float ScanlineTime (Time) ; 

ScanlineTime : units = "hours" ; 
short BlackBodylTemperature (Time, Number Of Channels ) ; 

BlackBodylTemperature : units = "degree_Celsius " ; 

BlackBodylTemperature: scale_f actor = 0 . 0099999998f ; 
short BlackBody2Temperature (Time, NumberOf Channels ) ; 

BlackBody2Temperature : units = " degree_Celsius " ; 

BlackBody2Temperature : scale_f actor = 0 . 0099999998f ; 
short Amplif ierGain (Time, NumberOf Channels ) ; 

Amplif ierGain : scale_f actor = O.OOlf ; 
short BlackBodylCounts (Time, NumberOf Channels ) ; 

short BlackBody2Counts (Time, NumberO f Channel s ) ; 

float CalibrationSlope (Time, NumberOf Channels ) ; 

CalibrationSlope : units = "watts /meter2 / steradian/micron" ; 
float Calibrationlntercept (Time, NumberOf Channels ) ; 

Calibrationlntercept : units = 

11 watts /me ter 2 / steradian/micron" ; 

float PixelLatitude (Time, AnchorlndexSize) ; 

PixelLatitude : units = " degree_north" ; 
float PixelLongitude (Time, AnchorlndexSize) ; 

PixelLongitude : units = " degree_east " ; 

float SensorZenithAngle (Time, AnchorlndexSize) ; 

SensorZenithAngle: units = "degree" ; 
float SensorAzimuthAngle (Time, AnchorlndexSize) ; 

SensorAzimuthAngle : units = "degree" ; 
float SolarZenithAngle (Time, AnchorlndexSize) ; 

SolarZenithAngle : units = "degree" ; 
float SolarAzimuthAngle (Time, AnchorlndexSize) ; 

SolarAzimuthAngle : units = "degree" ; 
float Aircraf tLatitude (Time ) ; 

Aircraf that itude : units = "degree__north" ; 
float Aircraf tLongitude (Time ) ; 

Aircraf tLongitude : units = " degree_east " ; 

float Aircraf tHeading (Time) ; 

Aircraf tHeading : units = "degree_true" ; 
float Aircraf tAltitude (Time) ; 

Aircraf tAltitude : units = "meters" ; 
float Aircraf tPitch (Time) ; 

Aircraf tPitch : units = "degree" ; 
short CalibratedData (Time , NumberOf Channels , NumberOf Pixels ) ; 

CalibratedData : units = "watts /meter2 /steradian/micron" ; 

CalibratedData : scale„f actor = O.f, O.lf, O.lf, 

0 . 0099999998f , 0 . 0099999998 f , 0 . 0099999998f , O.OOlf, 0 . 0099999998f , 

0 . 0099999998f , 0 . 0099999998f , 0 . 0099999998f , 0 . 0099999998f ; 

// global attributes: 

:title = "MODIS Airborne Simulator (MAS) Level-lB Data" ; 

: CreationDate = "27-Jan-94 10:39:54" ; 

: ExperimentName = "ASTEX " ; 

:FlightDate = "17 June 1992 " ; 

: FlightNumber = "92-107 " ; 

: ClockUsedToProcess = "INS Clock " ; 

: Flight LineNumber = 8 ; 

: TotalFlightLines = 13 ; 

: Credits = "Liam Gumley/design Paul Hubanks /modify (MODIS 

SDST ) " ; 
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: MASDataUsersGuideSource = "Paul Hubanks (MODIS SDST) 301 982 
3724 or hubanks@ltp.gsfc.nasa.gov" ; 

: ClockProblemDescription - "The MAS and INS (navigation) clocks 
often displayed different times. For the FIRE experiment a single clock 
offset was computed and and stored in the variable FIREClockOff set ( INS- 
MAS) Hours. The clock used in the FIRE experiment was the MAS Clock, and is 
stored in the array GreenwichMeanTime . For all post-FIRE experiments, 
either the MAS or INS clock may have been used. To determine this check 
the entry ClockUsedToProcess . This may be set to either MAS_Clock or 
INS__Clock. During a flight, for both FIRE and post-FIRE flights, time is 
computed from the recorded start time of each flight track, incrementing 
.16 seconds per scanline. This time is stored in ScanlineTime . Finally, 
for all post-FIRE experiements the entire string of INS and MAS clocktimes 
are stored in the arrays InsTime and MasTime" ; 

: FIREClockOff set (INS-MAS) Hours = -99. f ; 

: VIS/NIRCountSensitivityAdjustment = "Yes (See DataSetHeader ) 

: Sof twareVersion = "Version 3.0 " ; 

: CalibrationVersion = "ASTEX King 1.0 " ; 

: data_set = " MAS ASTEX ” ,* 

: data ^product = "straight-line flight tracks " ; 

:geog_flag = "c" ; 

: day_night__f lag = "d" ; 

: granule__version = 3 ; 

: metadata_version = " 3a" ; 

: producer_granule_id = "Version 3.0 ASTEX King 1.0 " ; 

: data_quali ty = "Fair: some noise in IR channels " ; 

: granule_size = 1 . 0884502e+08f ; 

:begin_date = "19920617 122120" ; 

:end_date = "19920617 123544" ; 

: lat_LL = 36.782097f ; 

: lon_LL = -24.800873f ; 

: lat_UL = 35 . 914349f ; 

: lon_UL = -23.128139f ; 

: lat_UR = 35.639751f ; 

: 1 on_UR = -23.356495f ; 

:lat_LR = 3 6. 501247 f ; 

: 1 on_LR = -25.026567f ; 

} 

7.2 Data Structure Term Definitions 
Dimensions: 

This section defines the dimensions and gives a brief description of each 
dimension contained in the data set. Dimensions are named integers which 
specify the shape of one or more multi-dimensional variables in the HDF file. 

Time = UNLIMITED ; // (5367 currently) 

The 'length' dimension of the data set. This dimension corresponds to elapsed 
time, and also corresponds with the incrementing of the scan line counter. 

NumberOf Channels = 12 ; 

The number of data channels. 
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NumberOfPixels = 716 ; 

The number of pixels (IFOVs) across a MAS scan. 

HeaderLength = 6160 ; 

The length of the ASCII data set header. The header is composed of blocks of 80 
characters, which contain no carriage return, line feed or similar record 
separators. The header contains a copy of the MAS instrument configuration file 
that was used in processing. 

AnchorlndexSize = 73 ; 

The number of geolocation anchor points per scan line. Each MAS scan line has 
geolocation data for every 10th pixel (to save space). The geolocation data is 

defined for pixel numbers 1, 10, 20, 30, 40, 680, 690, 700, 710 and 716 which 

makes a total of 73 geolocation anchor points. It should be noted that pixels 1 to 
358 are on the starboard (right) side of the aircraft, while pixels 359 to 716 are on 
the port (left) side of the aircraft. 

Variables: 

This section defines the type and size of the variables, along with a brief 
description of each, in the output data set. A variable represents a multi- 
dimensional array of values of the same type. The type definitions are from the 
C language, however the equivalent FORTRAN types are 


byte: 

BYTE 

8 bit data 

char: 

CHARACTER 

synonymous with byte 

short: 

INTEGERS 

16 bit integers 

long: 

INTEGERS 

32 bit integers 

float: 

REALM 

32 bit IEEE floating point 

double: 

DOUBLE PRECISION 

64 bit IEEE floating point 


Note that the listing of variable names which follows show the dimensions in C 
order. In FORTRAN, the order of dimensions is reversed. For example, the 
variable declared as CalibratedData(Time, NumberOfChannels, NumberOfPixels) 
in C would be declared as CalibratedData(NumberOfPixels, NumberOfChannels, 
Time) in FORTRAN. 

float Left50%ResponseWavelength(NumberOfChannels) ; 
Left50%ResponseWavelength:units = "microns" ; 

The wavelength corresponding to the left hand 50% response point of the 
spectral response function for each channel. 

float CentrallOO%ResponseWavelength(NumberOfChannels) ; 
CentrallOO%ResponseWavelength:units = "microns" ; 

The wavelength corresponding to the central 100% response point of the spectral 
response function for each channel. 

float Right50%ResponseWavelength(NumberOfChannels) ; 
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Right50%ResponseWavelength:units = "microns" ; 

The wavelength corresponding to the right hand 50% response point of the 
spectral response function for each channel. 

float SolarSpectrallrradiance(NumberOfChannels) ; 

SolarSpectralIrradiance:units = "watts/meter2/micron" ; 

SolarSpectrallrradiancertitle = "Sensor Weighted and Orbit Corrected Solar 
Spectral Irradiance" ; 

The sensor weighted and orbit corrected solar spectral irradiance at the top of the 
atmosphere for each channel. This is computed using the LOWTRAN7 database 
of solar spectral irradiance and the spectral response function for each channel by 

S W (X) = S S(X) . R(X) dX (12) 

J R(X) dX 


where 


S W (X) = sensor weighted orbit corrected exoatmospheric solar spectral 
irradiance, 

S(X) - orbit corrected exoatmospheric solar spectral irradiance, 

R(X) = sensor spectral response. 

short AnchorPointlndex(AnchorlndexSize) ; 

The pixel numbers for which geolocation information is defined. The 

geolocation data is defined for pixel numbers 1, 10, 20, 30, 40, 680, 690, 700, 710 

and 716 which makes a total of 73 geolocation anchor points. 

char DataSetHeader(HeaderLength) ; 

The data set header text. The header is composed of blocks of 80 characters, 
which contain no carriage return, line feed or similar record separators. The 
header contains a copy of the MAS instrument configuration file that was used 
in processing. This indicates, among other things, which spectral bands were 
enabled, which channels had 8 or 10 bit data, which channels were calibrated 
using the MAS blackbodies (IR channels), and slopes/intercepts for the 
visible /near-IR channels. 

short DataFrameStatus(Time) ; 

(LEVEL-0 MAS engineering data) 

Data Frame Status. 

Zero indicates good data. Bits are set to 1 to indicate errors as follows: 


Bit Position 

Meaning 

1 

Bad digit in scanline count 

2 

Bad digit in time code 

4 

Bad digit in thumbwheel data 

8 

Bad digit in reference temperature 

16 

Incorrect channel 

64 

Bad sync word 

128 

No match for end-of-frame code (last word) 
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sign Set to 1 if any other bit set 

short RunNumber(Time) ; 

(LEVEL-0 MAS engineering data) 

Not used. Was intended for operator to encode flight line number during flight. 

long ScanLineCounter(Time) ; 

(LEVEL-0 MAS engineering data) 

Scan line count, increments by 1 for each consecutive scan line. 

long Thumb Wheels witches (Time) ; 

(LEVEL-0 MAS engineering data) 

Data system thumbwheel switch settings formatted as YYFFFDDD, where YY is 
financial year, FFF is flight number, DDD is day of year. 

short ScanRate(Time) ; 

ScanRate:units = "hertz" ; 

ScanRate:scale_factor = O.lf ; 

(LEVEL-0 MAS engineering data) 

Scan rate in scans per second (x 10, nearest integer). 

long GreenwichMeanTime(Time) ; 

(LEVEL-0 MAS /INS engineering data) 

Greenwich Mean Time is either the MAS or INS time, selected at the time of 
processing. It is used to determine the time at the start of the straight line flight 
tracks and is formatted as HHMMSSS, where Hli is hours, MM is minutes, SSS 
is seconds and tenths of a second. The tenths of a second field is not currently 
recorded, therefore it is always zero. 

long MasTime(Time) ; 

(LEVEL-0 MAS engineering data) 

Time recorded by the MAS clock and is formatted as HHMMSSS, where HH is 
hours, MM is minutes, SSS is seconds and tenths of a second. The tenths of a 
second field is not currently recorded, therefore it is always zero. 

long InsTime(Time) ; 

(LEVEL-0 INS engineering data) 

Time recorded by the INS clock formatted as HHMMSSS, where HH is hours, 
MM is minutes, SSS is seconds and tenths of a second. The tenths of a second 
field is not currently recorded, therefore it is always zero. 

short S-Bendlndicator(Time) ; 

(LEVEL-0 MAS engineering data) 

S-bend indicator where 0 = no S-bend, 1 = S-bend. 

Not used. Was intended to indicate whether hardware geometric correction was 
enabled. 

short AircraftRollCount(Time) ; 

(LEVEL-0 MAS engineering data) 
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Aircraft roll count. Signed integer where positive is right, 0.03 degrees per count. 

long YearMonthDay(Time) ; 

Date formatted as YYYYMMDD where YYYY is year, MM is month, DD is day. 

float ScanlineTime(Time) ; 

ScanlineTime:units = "hours” ; 

Time of the current scanline in fractional hours, computed from the scanline 
number and scan rate (as described in the Section 5). ScanlineTime is accurate to 
the nearest scanline (0.16 seconds) whereas GreenwichMeanTime is accurate 
only to the nearest whole second. 

short BlackBodylTemperature(Time, NumberOfChannels) ; 
BlackBodylTemperature:units = "degree_Celsius" ; 
BlackBodylTemperature:scale_factor = 0.0099999998f ; 

(LEVEL-0 MAS engineering data) 

Black Body 1 (cool) thermal reference temperature (degrees C x 100). 

short BlackBody2Temperature(Time, NumberOfChannels) ; 
BlackBody2Temperature:units = "degree_Celsius" ; 
BlackBody2Temperature:scale_f actor = 0.0099999998f ; 

(LEVEL-0 MAS engineering data) 

Black Body 2 (warm) thermal reference temperature (degrees C x 100). 

short AmplifierGainfTime, NumberOfChannels) ; 

AmplifierGain:scale_factor = O.OOlf ; 

Gain setting for each channel (x 1000). 

short BlackBodylCounts(Time, NumberOfChannels) ; 

(LEVEL-0 MAS engineering data) 

Black Body 1 (cool) radiance digital counts. 

short BlackBody2Counts(Time, NumberOfChannels) ; 

(LEVEL-0 MAS engineering data) 

Black Body 2 (warm) radiance digital counts. 

float CalibrationSlope(Time, NumberOfChannels) ; 

CalibrationSlope:units = "watts/meter2/steradian/micron" ; 

Linear calibration slope for each channel for digital scene count to radiance 
conversion. 

float CalibrationIntercept(Time, NumberOfChannels) ; 

CalibrationIntercept:units = "watts/meter2/steradian/micron" ; 

Linear calibration intercept for each channel for digital scene count to radiance 
conversion. 

float PixelLatitude(Time, AnchorlndexSize) ; 

PixelLatituderunits = "degree_north" ; 
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Latitudes for pixels at geolocation anchor points. Latitude ranges from -90 
degrees at the South Pole to +90 degrees at the North Pole. 

float PixelLongitude(Time, AnchorlndexSize) ; 

PixelLongitude:units = "degree_east" ; 

Longitudes for pixels at geolocation anchor points. Longitude is zero at the 
Greenwich Meridian, and ranges from -180 degrees (West) to +180 degrees (East). 

float SensorZenithAngle(Time, AnchorlndexSize) ; 

SensorZenithAngle:units = "degree" ; 

Sensor zenith angle for pixels at geolocation anchor points. Defined as the 
zenith angle (degrees) of a vector from the sensor to the pixel (nadir sensor 
zenith angle = 0). 

float SensorAzimuthAngle(Time, AnchorlndexSize) ; 

SensorAzimuthAngle:units = "degree" ; 

Sensor azimuth angle for pixels at geolocation anchor points. Defined as the 
azimuth angle (degrees) clockwise from North of a vector from the pixel to the 
sensor. 

float SolarZenithAngle(Time, AnchorlndexSize) ; 

SolarZenithAngle:units = "degree" ; 

Solar zenith angle for pixels at geolocation anchor points. Defined as the zenith 
angle (degrees) of a vector from the pixel to the Sun. 

float SolarAzimuthAngle(Time, AnchorlndexSize) ; 

SolarAzimuthAngleiunits = "degree" ; 

Solar azimuth angle for pixels at geolocation anchor points. Defined as the 
azimuth angle (degrees) clockwise from North of a vector from the pixel to the 
Sun. 

float AircraftLatitude(Time) ; 

AircraftLatitude:units = "degree_north" ; 

Aircraft subpoint latitude (degrees, derived from linear regression of INS data in 
the flight track). 

float AircraftLongitude(Time) ; 

AircraftLongitude:units = "degree_east" ; 

Aircraft subpoint longitude (degrees, derived from linear regression of INS data 
in the flight track). 

float AircraftHeading(Time) ; 

AircraftHeadingnmits = "degree_true" ; 

Aircraft heading (degrees, derived from linear regression of INS data in the flight 
track). 

float Aircraft Altitude(Time) ; 

Aircraft Altitudeiunits = "meters" ; 
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Aircraft altitude (meters, derived from linear regression of INS data in the flight 
track). 

float AircraftPitch(Time) ; 

AircraftPitch:units = "degree" ; 

Aircraft pitch angle (degrees, derived from linear regression of INS data in the 
flight track). 

short CalibratedData(Time, Numb erOf Channels, Numb erOf Pixels) ; 
CalibratedDatarunits = "watts/meter2/steradian/micron" ; 
CalibratedData:scale_factor = O.f, O.lf, O.lf, 0.0099999998f, 0.0099999998f, 
0.0099999998f, O.OOlf, 0.0099999998f, 0.0099999998f, 0.0099999998f, 0.0099999998f, 
0.0099999998f ; 

Calibrated radiances for all pixels in all channels. 32 bit floating point radiance 
values are scaled to 16 bit integers to save space. Variable scaling factors are used 
for each channel. To retrieve real radiances, multiply the integer radiance value 
by the appropriate scaling factor for that channel. 


Global Attributes: 


This section describes the global attributes of the output data set. Global 
attributes contain descriptive information pertaining to the entire data set. 

: title = "MODIS Airborne Simulator (MAS) Level-IB Data" ; 

: CreationDate = "27-Jan-94 10:39:54" ; 

Date that this file was created (processed). 

: ExperimentName = "ASTEX " ; 

Name of the MAS mission. 

: FlightDate = "17 June 1992 " ; 

Date of the MAS flight. 

: FlightN umber = "92-107" ; 

Flight number assigned. 

: ClockUsedT oProcess = "INS Clock " ; 

Clock used to determine the start time of the straight line flight tracks. 

: FlightLineNumber = 8 ; 

Flight line number. 

:TotalFlightLines = 13 ; 

Total number of flight lines (straight-line flight tracks) during this day (flight). 
:Credits = "Liam Gumley/design Paul Hubanks/modify (MODIS SDST)" ; 
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: MASDataUsersGuideSource = "Paul Hubanks (MODIS SDST) 301-982-3724 or 

hubanks@ltp.gsfc.nasa.gov” ; 

: ClockProblemDescription = (Description of time keeping problems onboard ER- 

2 ); 

: FIREClockOffset (INS-MAS) Hours = -99 .f ; 

Offset between INS and MAS clocks for pre-Exabyte recording system LEVEL-0 
data only. Before the Exabyte recording system was brought into use in February 
1992 the MAS and INS data streams were recorded separately, each using it's own 
internal clock. During the FIRE mission, the INS clock was found to be leading 
the MAS clock by anywhere from 30 to 90 seconds. This offset had to be corrected 
to ensure accurate geolocation, so the INS clock time was adjusted in post- 
processing to match the MAS clock time. This variable records the time 
subtracted from the INS clock time. After February 1992 the Exabyte LEVEL-0 
recording system encoded the MAS and INS in one synchronous data stream, so 
this problem was eliminated. 

: VIS/NIRCountSensitivityAdjustment = "Yes (See Data Set Header) ” ; 

Flag to indicate form of adjustment (if any) to VIS/NIR digital counts to 
compensate for instrument temperature dependent changes in sensitivity. See 
the section on calibration for details. NONE indicates no adjustment. If used, 
the form of the adjustment is listed in the instrument configuration file which is 
stored in DataSetHeader. 

: SoftwareVersion = "Version 3.0 " ; 

Version of the software used to create this HDF file (MAS LEVEL-1B data 
granule). 

: CalibrationVersion = "ASTEX King 1.0 " ; 

Version of the calibration used in the processing of this particular HDF file (MAS 
LEVEL-1 B data granule). 

: data_set = "MAS ASTEX " ; 

Data set name (DA AC required metadata). 

: data_product = "straight-line flight tracks " ; 

Data product name (DA AC required metadata). 

: geog_flag = "c" ; 

geographic flag (DA AC required metadata) c=constrained. 

: day_night_flag = "d" ; 

Day /Night flag (DA AC required metadata). 

: granule_version = 3 ; 

Version of the granule. Shows the number of time of times this data has been 
processed in the past (DA AC required metadata). 
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: metadata_version = " 3a" ; 

Version of the metadata included in this particular HDF file (MAS LEVEL-1B 
data granule). Allows indication of an update to the metadata when no update is 
done to the calibrated radiances. 

: producer_granule_id = "Version 3.0 ASTEX King 1.0 " ; 

A unique ID assigned by the data producer (DA AC required metadata). 

: data_quality = "Fair: some noise in IR channels " ; 

Data quality assessment (DA AC required metadata). 

: granule_size = 1.0884502e+08f ; 

Approximate size of the granule in bytes (DAAC required metadata). 

: begin_date = "19920617 122120" ; 

Date and time of the start of the flight track (granule) (DAAC required metadata). 

: end_date = "19920617 123544" ; 

Date and time of the start of the flight track (granule) (DAAC required metadata). 

rlatJLL = 36.782097f ; 

:lon_LL = -24.800873f ; 

:lat_UL = 35.914349f ; 

:lon_UL = -23.128139f ; 

:lat_UR = 35.639751f ; 

:lon_UR = -23.356495f ; 

:lat_LR = 36.501247f ; 

:lon_LR = -25.026567f ; 

Latitude and longitude of the corner points of the granule (DAAC required 
metadata). 

} 
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8. Obtaining Summary and Browse Information 


MAS flight summary information and browse data are available via anonymous 
FTP at ltpiris2.gsfc.nasa.gov in directory 

pub/MAS/data 

This site may be accessed as shown: 

% ftp ltpiris2.gsfc.nasa.gov 
Name : anonymous 
Password: guest 
ftp> cd pub/MAS/data 
ftp> dir 


-rw-r--r— 

1 

11421 

20 

868 

Oct 

14 

09 : 16 

00 readme 

drwxr-xr-x 

13 

11421 

20 

512 

Sep 

1 

1993 

ASTEX 

drwxr-xr-x 

3 

11421 

20 

512 

Mar 

21 

10:03 

BOREAS -hall 

drwxr-xr-x 

2 

11421 

20 

512 

Oct 

5 

1993 

CEPEX 

drwxr-xr-x 

18 

11421 

20 

512 

Oct 

14 

09:25 

FIRE 

drwxr-xr-x 

11 

11421 

20 

512 

Oct 

20 

10:14 

SCAR-A 

drwxr-xr-x 

3 

11421 

20 

512 

Jul 

9 

1993 

TOGA 


The readme file contains information on MAS LEVEL-1B data processing status, 
and data availability for each experiment. In addition, directories have been 
created to hold summary and browse information for each MAS experiment. 
For example, to obtain browse and summary data for a particular ASTEX flight, 
enter the ASTEX directory by typing: 


ftp> cd ASTEX 
ftp> dir 


-rwxr-xr-x 

1 

11421 

20 

403 

Aug 

12 

1993 

0 0 README . DOC 

drwxr-xr-x 

2 

11421 

20 

512 

Apr 

5 

12:50 

02 jun92 

drwxr-xr-x 

2 

11421 

20 

1536 

Apr 

5 

12:55 

04 j un92 

drwxr-xr-x 

2 

11421 

20 

1536 

Apr 

5 

12 : 58 

08 j un92 

drwxr-xr-x 

2 

11421 

20 

1536 

Feb 

2 

09:17 

17 j un92 

drwxr-xr-x 

2 

11421 

20 

1536 

Oct 

18 

14 : 58 

18 jun92 

drwxr-xr-x 

2 

11421 

20 

1536 

Oct 

18 

15:01 

19 jun92 

drwxr-xr-x 

2 

11421 

20 

1024 

Feb 

16 

09:42 

21 jun92 

drwxr-xr-x 

2 

11421 

20 

1024 

Feb 

16 

13:14 

22 jun92 

drwxr-xr-x 

2 

11421 

20 

1024 

Oct 

18 

15 : 06 

23 jun92 

drwxr-xr-x 

2 

11421 

20 

1536 

Feb 

9 

06 : 55 

2 9 may 9 2 

drwxr-xr-x 

2 

11421 

20 

1024 

Oct 

18 

15 : 09 

3 lmay 9 2 


Within each experiment directory there are a number of subdirectories. Each 
subdirectory contains data recorded on a particular day. The sub-directories are 
named with the day, month, and year to enable easy selection of a particular day 
of interest. Each sub-directory contains the flight summary (*.sum) and browse 
(*.gif) data for one MAS flight during the experiment. For example, enter the 
directory for 17 June 1992 by typing: 


% cd 17 jun92 
% dir 

-rw-r— r— 1 11421 20 2307 Feb 2 09:12 92169. sum 

-rw-r — r — 1 11421 20 39022 Feb 2 09:11 9216901i.gif 
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-rw-r- 

-r- - 

1 

11421 

20 

41241 

Feb 

2 

09 : 11 

9216901v.gif 

-rw-r- 

-r- - 

1 

11421 

20 

75311 

Feb 

2 

09:12 

9216902i.gif 

-rw-r- 

-r-- 

1 

11421 

20 

64532 

Feb 

2 

09:12 

9216902v. gif 

-rw-r- 

-r-- 

1 

11421 

20 

18146 

Feb 

2 

09 : 12 

9216903i.gif 

-rw-r- 

-r-- 

1 

11421 

20 

19379 

Feb 

2 

09:12 

9216903v.gif 

-rw-r- 

-r-- 

1 

11421 

20 

49289 

Feb 

2 

09:12 

9216904i.gif 

-rw-r- 

-r-- 

1 

11421 

20 

43389 

Feb 

2 

09 : 12 

9216904v.gif 

-rw-r- 

-r-~ 

1 

11421 

20 

98800 

Feb 

2 

09:12 

9216905i.gif 

-rw-r- 

-r-- 

1 

11421 

20 

98075 

Feb 

2 

09 : 12 

9216905v.gif 

-rw-r- 

-r-- 

1 

11421 

20 

41049 

Feb 

2 

09 : 12 

9216906i.gif 

-rw-r- 

-r-- 

1 

11421 

20 

28162 

Feb 

2 

09 : 12 

9216906v. gif 

-rw-r- 

-r-- 

1 

11421 

20 

18944 

Feb 

2 

09:12 

9216907i.gif 

-rw-r- 

-r — 

1 

11421 

20 

18570 

Feb 

2 

09 : 12 

9216907v.gif 

-rw-r- 

-r — 

1 

11421 

20 

198545 

Feb 

2 

09:12 

9216908i.gif 

-rw-r- 

-r — 

1 

11421 

20 

182063 

Feb 

2 

09 : 12 

9216908v.gif 

-rw-r- 

-r-- 

1 

11421 

20 

173019 

Feb 

2 

09:12 

9216909i.gif 

-rw-r- 

-r- - 

1 

11421 

20 

162709 

Feb 

2 

09 : 12 

9216909v.gif 

-rw-r- 

-r — 

1 

11421 

20 

115572 

Feb 

2 

09 : 12 

9216910i.gif 

-rw-r- 

-r- - 

1 

11421 

20 

111105 

Feb 

2 

09:12 

9216910v.gif 

-rw-r- 

-r-- 

1 

11421 

20 

123971 

Feb 

2 

09:12 

9216911i.gif 

-rw-r- 

-r-- 

1 

11421 

20 

123051 

Feb 

2 

09:12 

9216911v.gif 

-rw-r- 

-r-- 

1 

11421 

20 

137309 

Feb 

2 

09 : 12 

9216912i.gif 

-rw-r- 

-r-- 

1 

11421 

20 

136656 

Feb 

2 

09 : 12 

9216912v . gif 

-rw-r- 

-r-- 

1 

11421 

20 

184931 

Feb 

2 

09 : 12 

9216913i.gif 

-rw-r- 

-r-- 

1 

11421 

20 

161166 

Feb 

2 

09 : 12 

9216913v . gif 

The flight summary file (92169. sum) is an 

ASCII 

text summary (and therefore 

must be transferred in 

ASCII mode 

in FTP) of the 

flight lines processed for that 


flight. The name is formatted as yyddd.sum where yy is the year and ddd is the 
day of year. The contents of this file are shown below: 


MODIS AIRBORNE SIMULATOR (MAS) FLIGHT LINE INFORMATION FOR 17-JUN-1992 




START 

OF FLIGHT 

LINE 



LINE 

TIME 

LAT 

LON 

SOLAR 

START 


HH :MM: SS 

DEG 

DEG 

ZEN 

AZIM 

HEADING 

1 

11:27:39 

37.214 

-24.754 

31.5 

106.6 

126.35 

2 

11:32:26 

36.831 

-24.421 

30.2 

107.5 

198.62 

3 

11:45:29 

36.155 

-24.434 

27 . 5 

109.8 

126.42 

4 

11:48:26 

36.031 

-24.079 

26.6 

110.8 

96.75 

5 

11:52:30 

35.968 

-23.536 

25.4 

112.5 

107.87 

6 

12:05:43 

35.630 

-22 . 662 

22.2 

117.5 

296.53 

7 

12:11:41 

35.957 

-22.895 

21.5 

120.2 

116.74 

8 

12:21:21 

35.777 

-23.243 

20.0 

123.3 

304.01 

9 

12:41:04 

36.851 

-25.025 

18 . 6 

132.1 

125.78 

10 

12:55:40 

35.802 

-23.544 

15.0 

142.9 

192.26 

11 

13:06:34 

35.058 

-24.057 

13.4 

148.3 

20.87 

12 

13:14:28 

35.886 

-23 .716 

13.3 

158.4 

17.88 

13 

13:25:52 

37.101 

-23.440 

13.8 

171.6 

296.18 


END OF FLIGHT LINE 


TIME 

LAT 

LON 

SOLAR 

SCAN 

HH:MM: SS 

DEG 

DEG 

ZEN 

AZIM 

LINES 

11:30:35 

37.016 

-24.435 

30.6 

107.3 

1098 

11:38:12 

36.213 

-24.686 

29.1 

107.6 

2169 

11:46:48 

36.065 

-24.292 

27 .1 

110.1 

495 

11:52:05 

35.977 

-23.590 

25.5 

112.4 

1370 

11:59:06 

35.747 

-22.707 

23.5 

115.3 

2473 

12:08:16 

35.754 

-22.996 

22.0 

118.3 

963 

12:12:56 

35.890 

-22 .750 

21.1 

120.8 

464 

12:35:39 

36.642 

-24.914 

19.2 

129.0 

5367 

12:54:08 

35.952 

-23.595 

15.4 

141.9 

4898 

13:02:50 

34.996 

-23.779 

13.6 

146.1 

2693 

13:14:13 

35.865 

-23.724 

13.3 

158.1 

2872 

13:24:06 

36.933 

-23.359 

13.7 

170.2 

3619 

13:37:23 

37.695 

-24.975 

14.3 

177.0 

4325 


NUMBER OF FILES FOR THIS FLIGHT = 13 
TOTAL NUMBER OF SCAN LINES = 32806 
DATE THESE FILES WERE PROCESSED = 25-Jan-94 
DATE THIS LIST WAS CREATED = 25-Jan-94 


The *.gif files shown in the previous directory listing are the browse images in 
Graphics Interchange Format (GIF) for each flight line. These images contain 
every 4th pixel on every 4th line from one short-wave visible and one long- 
wave infrared channel. They are contrast stretched and converted to GIF for 
distribution, and as such do not contain quantitative radiance information. 
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There are 2 GIF image files per flight line, named 92170??v.gif and 9217077i.gif, 
where ?? is the flight line number, v denotes visible (VIS) imagery, and i denotes 
infrared (IR) imagery. The date corresponds to year and day of year as mentioned 
previously. The VIS and IR channels are selected such that the VIS channel is at 
around 0.68 microns, and the IR channel is at around 12 microns. 

Example browse images from the MAS flight on 17 June 1992 are shown on the 
following page (Figure 1). These were acquired over the Atlantic Ocean, off the 
coast of the Azores Islands. Stratiform cloud features are visible, with marine 
stratocumulus in the upper region of the image and continental stratocumulus 
in the lower region of the image. A cold cirrus cloud formation can be seen (as a 
black feature) in the infrared browse image in the upper portion. This feature is 
not distinct in the visible browse image. 

Public domain or shareware software packages to view GIF imagery on various 
systems are available from the following anonymous FTP sites: 

IBM/PC 

oak.oakland.edu, pub/msdos/ gif/ cshwl01a.zip 

Macintosh 

mac.archive.umich.edu, mac/graphics/ graphicsutil/ gifconverter2.37.cpt.hqx 
Unix/X 

ftp.x.org, contrib /xloadimage.4.1 .tar.gz 
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Visible (0.68 \im) browse image 


Infrared (11.95 jam) browse image 


Figure 1. MAS Browse Images, Flight Line 11, 17 June 1992 
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9. MAS Contact List 


Below is a list of persons involved with the processing and analysis of MAS data. 

Paul Hubanks, MAS LEVEL-1B Processing Manager 

Research and Data Systems Corporation 

7855 Walker Drive, Suite 460 

Greenbelt MD 20770, USA 

(301) 982-3724, phubanks@ltpsun.gsfc.nasa.gov 

Liam Gumley, MAS LEVEL-1 Processing System Designer 
Code 913, NASA Goddard Space Flight Center 
Greenbelt MD 20771, USA 
(301) 286-8789, gumley@climate.gsfc.nasa.gov 

Tom Arnold, MAS VIS/NIR Calibration 
Code 913, NASA Goddard Space Flight Center 
Greenbelt MD 20771, USA 
(301) 286-4805, arnold@climate.gsfc.nasa.gov 

Michael King, MAS Principal Investigator 
Code 900, NASA Goddard Space Flight Center 
Greenbelt MD 20771, USA 
(301) 286-8228, king@climate.gsfc.nasa.gov 

Paul Menzel, MAS Principal Investigator 
N O A A / NESDIS 
1225 W. Dayton St. 

Madison WI 53706 USA 

(608) 263-4930, paulm@ssecmail.ssec.wisc.edu 

Chris Moeller, MAS Investigator; calibration, geolocation, MCIDAS 
Space Science and Engineering Center, UW-Madison 
1225 W. Dayton St. 

Madison WI 53706, USA 

(608) 263-9597, chrism@ssecmail.ssec.wisc.edu 

Jeff Myers, MAS LEVEL-0 Data Distribution 

High Altitude Missions Branch, NASA Ames Research Center 

Mail Stop 240-6 

Moffett Field CA 94035, USA 

(415) 604-6252, jmyers@msmail.arc.nasa.gov 

Pat Grant, MAS Engineering 

High Altitude Missions Branch, NASA Ames Research Center 

Mail Stop 240-6 

Moffett Field CA 94035, USA 

(415) 604-6393, pgrant@msmail.arc.nasa.gov 
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Appendix A. MAS Spectral Bands 
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Gaussian curves fitted to estimated 50% response points. 
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